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Abstract 

We study the distribution of the color fields due to a static quark-antiquark 

pair in SU(2) lattice gauge theory. We find that the London penetration length 

measured after Abelian projection in the Abelian Covariant gauge (Maximal 

Abelian gauge) agrees with the one obtained without gauge fixing. Moreover 

the penetration length scales according to asymptotic freedom. We put out a 

simple relation between the penetration length and the string tension. 
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To understand the non perturbative phenomenon of color confinement in non Abelian 
gauge theories, G. 't Hooft [|TJ and S. Mandelstam proposed a model known as dual 
superconductor model. The physical grounds of the model lie in the theory of superconduc- 
tivity 0. In the 't Hooft's formulation the dual superconductor model is elaborated in the 
framework of the Abelian projection ||. After a particular gauge has been fixed, reducing 
the symmetry to that of the maximal Abelian (Cartan) subgroup, the non Abelian gauge 
theory is described in terms of Abelian projected gauge fields ("photons"). 

In this scenario there are also color magnetic monopoles whose condensation should 
cause the confinement of all particles which are color electrically charged with respect to the 
above photons. Of course this definition of the color magnetic monopoles does depend on 
the Abelian gauge fixing. However the physics, i.e. the monopole condensation, should be 
independent on the gauge fixing. 

In this paper we analyze the finger-print of the dual superconductor hypothesis, namely 
the dual Meissner effect. To do this we analyze the distribution of the color fields due to a 
static quark-antiquark pair in SU(2) lattice gauge theory. Following Ref. one can explore 
the field configurations produced by the quark-antiquark pair by measuring the connected 
correlation function: 

(tr(WLUpLt)) 1 (tr(Up)tr(W)) 
PW (tr(W)> 2 (tr(W)> ' { ' 

where the plaquette Up = U^ v {x) in the (p,v) plane is connected to the Wilson loop W 

through the Schwinger line L. Note that the correlation function ([I]) is sensitive to the field 

strength rather than to the square of the field strength ||: 



Pw — -* a 2 g 



(F^ qq -(F^) . (2) 



According to Eq.(H) we define the color field strength as: 

Fnu{x) = ~yPw(x) . (3) 

By varying the distance and the orientation of the plaquette U p with respect to the Wilson 
loop W, one can scan the color field distribution of the flux tube. 



We performed numerical simulations on 16 4 , 20 4 , and 24 4 lattices in the range 2.45 < 
P < 2.7. We used Wilson loop L w x L w , L w = L/2 — 2, L being the lattice size. In order 
to eliminate the uninteresting short range quantum fluctuations we cooled our statistical 
samples. In agreement with previous studies |5|U , the connected correlator pw turns out to 
be sizeable when the plaquette and the Wilson loop are parallel. Thus, the component of 
the chromoelectric field parallel to the line joining the static charges is sizeable, while the 
other components of the chromoelectric field and the chromomagnetic field are much smaller. 
Moreover the longitudinal chromoelectric field is almost constant along the quark- ant iquark 
line, and it decreases rapidly as the transverse distance is increased. 

Remarkably enough, in a previous study in the Abelian covariant gauge (maximal 
Abelian gauge) we found a similar behaviour. 

In the continuum the Abelian covariant gauge corresponds to impose the constraints 

D IM Afa)=0 (4) 

where = A], ± iA 2 „. and D„ is the A^ '-covariant derivative. On the lattice the con- 
straints (|]) are implemented by by maximizing || 

R = E [°Mx)a 3 Ul(x)] , (5) 

where U^(x) are the gauge- fixed links: 

U^x) =V(x)U^x)V\x + fi) . (6) 

It is straightforward to check that the residual gauge invariance is the U(l) group with 
transformations g(x) = exp [ia^9{x)\. 

In the dual superconductor scenario the long range properties of the gauge system are 
encoded into the Abelian fields. On the lattice the Abelian fields are defined through the 
Abelian projected links Ufi(x) 0: 

Ufa) = diag [e< (x) , e -<M] , 9 fa) = arg [^(x)] ^ . (7) 



In Ref. we considered the Abelian projected correlator pw- 

_ (tr(W A U A )) l (tr(U A )tr(W A )) 
Pw (tr(W A )) 2 (tr(W A )) ' {) 

We found that the longitudinal chromoelectric field in the middle of the flux tube can be 
fitted according to 

E,M = l±K (i) , * ± > (9) 

where xj_ is the transverse distance. Kq is the modified Bessel function of order zero, <3> is 
the external flux, and A is the London penetration length. 

Equation (H) is a straightforward consequence of the dual superconductor hypothesis. 
Indeed, the transverse shape of the longitudinal chromoelectric field is the dual version of 
the Abrikosov vortex field distribution. However, the chromoelectric field has been obtained 
through the gauge-dependent correlator p^- So the physical meaning of Eq. (|) is not clear. 
On the other hand, if the confinement is realized by the dual superconductor mechanism, 
then Eq.(^) should hold even if we use the gauge invariant correlator py/. To check this 
point, we measured Ei(x±_) in the middle of the flux tube both in the Abelian covariant 
gauge and without gauge fixing by varying xj_ up to transverse distance L w . 

Remarkably enough, we find that Eq. (Q) holds for both definition of the longitudinal 
chromoelectric field ||. We obtain a rather good fit to the data (x 2 /d.o.f < 1) if x± > 2 
(in lattice units); this is reasonable since Eq. (H) applies in the region A ^> £, where £ is the 
coherence length (more on this later on). 

In Figure 1 we report our results for the inverse of the penetration length p = 1/A. For 
the gauge invariant correlator pw we collected up to 100 configurations with 8 cooling steps 
(we, however, checked that the results are stable up to 12 cooling steps) at each (3. We 
performed measurements for 9 different values of (3 in the range 2.45 < (3 < 2.7 (for figure 
readability not all the values are displayed). 

In the Abelian covariant gauge we performed 100-500 measurements for 7 different values 
of (3 in the same range. The gauge is fixed iteratively via the overrelaxation algorithm of 
Ref. |l0l with the overrelaxation parameter uj = 1.7. 



Figure 1 suggests that both penetration lengths scale according to two-loop asymptotic 
scaling law. We find fi/ Aj^ = 8.9(3) (x 2 /d.o.f. ~ 2) for the gauge invariant correlator, and 
fi/AjYs — 8.3(7) (x 2 /d.o.f. ~ 0.4) in the Abelian covariant gauge. A striking consequence 
of Fig.l is that both definitions of the penetration length agree. The overall fit of the data 
gives (dashed lines) 

^ = 8.8(3) (10) 



with % 2 /d.o.f. ~ 1.8. Thus our results strongly suggests that the London penetration length 
is a physical gauge-invariant quantity. 

In Figure 2 we display the external flux $ versus (3. Naively we expect that 

$ = 1. (11) 

Indeed in the Abelian covariant gauge, we see that the data are quite close to Eq. ([LI]). 
On the other hand, without gauge fixing the external flux is strongly affected by lattice 
artefacts which seem to disappear by increasing (3. Thus we are led to suspect that the 
external flux gets renormalized by irrelevant operators, whose effects are strongly reduced 
in the Abelian covariant gauge. As we will argue in a moment, these effects are responsible 
for the non-scaling behaviour of the string tension. Indeed, the string tension can be defined 
as the energy stored into the flux tube per unit length. By observing that the longitudinal 
chromoelectric field is almost constant along the flux tube, we obtain from Eq. (H): 

_ $ 

a=-;=\i. (12) 
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Now, we have shown that \x almost scales according to asymptotic freedom. Whence the non- 
scaling behaviour of the string tension is due to the lattice renormalization of the external 
charge. However, Eq. (O) allows us to get rid of these effects by putting $ = 1 into Eq. (O) : 



~° = J k- (13) 
In Figure 3 we show Eq. (0). For comparison we display also the string tension extracted 
from the Wilson loops (full points). Due to low statistics and limited lattice size our ex- 
timation of the string tension is poor. However, it is gratifying to find that the linear 



extrapolation to the continuum limit of our data is compatible, albeit with large error, with 
the one obtained on larger lattices JIT) (full star): 

= 1.79 ±0.12. (14) 



The string tension defined by Eq. (0) (open points) seems to scale for f3 > 2.45 to a value 
which is consistent with (|i~4|). Indeed we find (dashed lines) 

= 1.76 ±0.15 , (15) 



where the quoted error include our extimation of the systematic error. 

We would like to comment on the validity of Eq. (|l"3D . In order to obtain Eq. (|T3|) we 
extrapolated Eq. (§) up to x± = 0. This extrapolation leads to a negligible error in fllTf ) if 
A/£ > 2. So it is important to have an extimation of £. The coherence length is determined 
by the monopole condensate, the order parameter for the confinement. Unfortunately we 
do not yet have at disposal a good definition of the order parameter. 

Recently two different groups [0,0] give an extimation of the coherence length by using 



monopole currents defined following DeGrand and Toussaint |TJ]]. However we feel that 
one can not rely on this approach. As a matter of fact, the authors of Ref. [|HJ argued 
convincingly that the DeGrand- Toussaint monopole density is not the order parameter for 
the confinement. 

Let us conclude by stressing the main results of the paper. We found that the London 
penetration length in the dual Meissner effect is a physical quantity. Moreover we put out 
a very simple relation between the string tension and the penetration length. A remarkable 
consequence of our Eq. (|j) is that the long range properties the SU(2) confining vacuum 
can be described by an effective Abelian theory. In addition, after fixing the gauge with 
the constraints (f|), it seems that the degrees of freedom which are not relevant to the 
confinement get suppressed. In this respect the situation looks like the Bardeen-Cooper- 
Schrieffer (BCS) theory of superconductivity. As it is well known |3|], in the BCS theory 
one considers a reduced Hamiltonian which, however, breaks the electromagnetic gauge 

6 



invariance. Nevertheless, the reduced BCS Hamiltonian, by retaining the degrees of freedom 
relevant to the superconductivity, gives the correct explanation of the Meissner effect. 

Let us conclude by stressing that that the most urgent problem to be addressed in the 
future studies is the reliable extimation of the coherence length [|1(| . 
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FIGURES 

FIG. 1. The inverse of the penetration length in units of Ajjg as a function of aA-j^. Unfixed 
gauge: circles L = 16, triangles L = 20, square L = 24. Abelian covariant gauge: crosses L = 16, 
diamond L = 20. 



FIG. 2. The external flux versus (3. Symbols as in Fig.l. 

FIG. 3. The square root of the string tension in units of A-^g versus aAjj-g. Open symbols 
refer to Eq. (jl3|). Full star is the extrapolated continuum limit. 
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